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ABSTRACT 

Useful relationships were formulated to describe variations of the diffuse spectral reflectance 
in terms of vegetation canopy variables, such as biomass. The relationships were based on 
the solution of the two-stream approximation of the radiative transfer equation. Out of 
the lengthy original e.xpression of the diffuse reflectance fonnula, simple working equations 
were derived by employing characteristic parameters, which are independent of the canopy 
coverage and identifiable by field observations. Tire typical asymptotic nature of reflectance 
data that is usually observed in biomass studies was clearly explained. The usefulness of 
the simplified equations was demonstrated by the exceptionally dose fit of the theoretical 
curves to two separately acq'rired data sets for alfalfa and shortgrass prairie canopies. 
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INTRODUCTION 


Agricultural remote sensing research in recent years has been largeiy concerned with develop- 
ing a fundamental quantitative understanding of the relationships between spectral responses 
and vegetation scene factors and how they are related to important agronomic parameters, 
such as plant biomass and crop yieid. Such quantitative relationships are essential for ex- 
tracting useful information from remotely sensed data for applications, such as forage 
management and pre-harvest prediction of crop yield. 

Numerous empirical relationships have been proposed for crop canopy assessments. Most 
of these remote sensing techniques rely on red and near-infrared reflectance or radiance 
ratios, which have been summarized by Tucker (1979). There is considerable evidence for 
a variety of cover types that red and photographic infrared spectral data are highly sensitive 
to the projected green leaf area index or green leaf biomass (Deering, 1978; Tucker, 1979; 
Holben et ah, 1980). Additionally, similar techniques have been found useful for indirectly 
assessing drought stress (Thompson and Wehmanen, 1979) and evapotranspiration (Wiegand 
et ah, 1979). 

Limitations and inconsistencies in the spectral relationships among the various cover types 
and conditions have prevailed, however, due to their lack of a theoretical foun'.lation. 

Tucker (1980) concluded that their utility in assessing standing crop biomass is tied to the 
relationship of the green leaf area index to the standing crop biomass for the cover type 
in question, and thus tliese relationships to standing crop biomass are not temporally 
consistent. Failure to take spectral measurements near solar noon results in additional 
inconsistencies in empirical models due to solar zenith angle effects (Duggin, 1977; 

Kriebel, 1978; Kimes et ah, 1980). 



Tills investigation was directed toward developing a basic relationship between vegetation 
canopy variables and dittuse spectral reflectance based on the radiative transfer theory. 

The Kubelka-Munk model (Allen and Richardson, 1968) was e.xtended to account for 
anisotropic diffusion of light within the canopy. Tlie goal of this study was to establish 
a practical procedure for analysis of biomass/reflectance data. 

MODIFICATION OF KUBELKA-MUNK MODEL 
Variations of the monochromatic diffuse radiation within vegetative canopies have been 
described by means of a two-parameter concept involving coefficients of absoiption and 
scattering (Allen and Richardson, 1968). This treatment is known as the Kubelka-Munk 
(KM) theory, which is applicable to homogeneous, perfectly diffusing media with light- 
absorbing and light-scattering elements. Such a treatment is seldom exact, but light intensity 
passing half-transparent materials has often been well approximated by the theory (Wend- 
landt and Hecht, 1966; Kortiim, 1969). 

Most vegetative canopies consist of several distinct components that result in anisotropic 
canopy reflectance. Examples are the reflectance differences between upper and lower 
surfaces of many plant leaves (Gates et al„ 1965) and bidirectional scattering effects of 
individual leaves (Breece and Holmes, 1971). Another example is inhomogeneous distribu- 
tions of leaf orientation. Tlie traditional two-parameter representation of the radiative 
transfer equation is certainly inadequate to take into account such phenomena. For this 
reason, the KM equations were extended by employing the two different sets of absorption 
and scattering coefficients as follows: 

dE- 

_ (a + y ) E. + 7.J. E+ (I) 

d (pz) 

dE4. 

d (pz) = - (a+ + 7 +) E+ + 7 _ E_ 


n 


( 2 ) 


where 


E_ = monochromatic llglit intensity in the downward direction 
= monochromatic liglit intensity in tlie upward direction 
z ° distance from tlie canopy top (negative in the downward direction) 
a_ s absorption coefficient associated with E_ 

= absorption coefficient associated with E+ 

= scattering coefficient associated with E_ 

74 . s scattering coefficient associated with 

p = density of plant canopy variable (e.g., biomass per unit volume). 

Here, d (pz) is the differential of biomass (biomass per unit area in a vertical distance 
segment dz). The variable “pz” can be converted into other canopy variables, such as 
leaf area index. The clioice of the appropriate parameter is often a matter of convenience 
in describing relatlonsliips between measured quantities. However, biomass is preferable to 
many otlier parameters, since the volume scattering by randomly oriented leaf elements is 
more appropriate for most plant canopies tlian multiple scattering in stratified leaf layers. The 
notations are simplified by omitting subscripts for the wavelength dependence of the 
variables and parameters. 

In this formulation the backsvard scattering of only the diffuse liglit is taken mi to account 
so that the problem remains one-dimensional (Fig. 1). It is a.ssumcd that there are negligible 
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Fig. 1. One-dimensional Plant Canopy Model 



contributions due to the emissivlty of any substances in the wavelengtlt considered. Otlter 
underlying assumptions are that the coefficients of bulk absorption and scattering are greater 
than zero. These assumptions are needed to avoid the problem of ambiguity, which occurs 
in cases where any coefncients vanish or become negative. The set of the present governing 
equations is a modification of tlie KM model, in wlilcli « a^. and 7 „ “ 7 ^. 

The gener.*'] solution set to the coupled differential equations is given by 


E_ = + A^e'^2^* (3) 

+ B2e">2^^ (4) 

where 

m, = D+K (5) 

m. « D-K (6) 

D == (a_ - + 7 .- “ 

K = (D^ + a_ + 7 ^ 7+ + (8) 


and Aj, A 2 , Bj and Bj are constants to be determined by boundary conditions. It can be 
sliown without difficulty tliat 

nij > 0 

m, < 0 

and also that, if a_ = and 7_ = 7 ^ , 

D =0 
m, = -m,. 

Convention shows in this type of problem that a common boundary condition is 

E_ = E^ at z = 0, (B.C. 1) 

where Eq is the intensity of the incident monochromatic light. 

Large Biomass Canopy Solution 

An important physical insiglit to the problem can be learned by examining the following 
hypothetical boundary condition: 
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finite at z « 


(B.C. *:) 


E_ and E^: 

vvlien the canopy heiglit or biomass is infinitely large. This hypotuetical case yields tlie 


solution set given by 

^ (9) 

E+ « RyE« (10) 

where 

Sj, = (a^ + a+ + T„ + 7+)/(2 7+) (12) 

“ K/7^ (13) 


and z < 0. In the present expression, Ry is the diffuse reflectance of the infinitely tall 
(or dense) canopy, which is a hypothetical situation. It is the reflectance of the vegetation 
canopy by which the influence of the background reflectance is eliminated. Tlic so-called 
KM relationship for diffuse reflectance (Park, 19S0); 

■V = > 7 - y<7>^ ^ 


a 

7 


(1 - Ry)- 
2 R„ 


(14b) 


is the special case of Eq. 11 where a = and y = 7 ^. In Eqs, II and 14a 

it is important to note that the diffuse reflectance is a function of ratio parameters, 
u^/7+> ^+/ 7 + and 7 _/ 7 ^ or a/ 7 . In Eq. 9, iHj is the mass attenuation coefficient of 
light, and the product nijP is the reciprocal of the penetration depth when the canopy 
is infinitely tall or dense. Eq. 9 is the typical form of the Beer-Lambert law (Wendlandt 
and Hecht, 1966). 


Model Solution 

Tlie background effects cannot be excluded in most canopy reflectance measurements. To 
account for the soil back, ground effects in the model of canopy height H, tlie second 
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boundary condition is substituted by 

E+ « RjjE. at z -II (B.C. 3) 

wljcre R{j is tlie reflectance of tlie background surface under Uu* canopy. Hie solution 
to Eqs, 1 and 2 satisfying Boundary Conditions 1 and 3 is given by 

- (Rv-Rb)'“= USa) 

Rn ■ ‘15b) 

" Rb'Rv'™''’‘“'^’’'-(Rv-Rb>— (1&) 

V Ri» 


I 

“f" I'^N + K„ - Rb> (Sn - 

- - Rn - Rb> '■‘■•n * (|<>m 

wliere 

C = (R^, - R^)»'"!PH (1711) 

N 

= (Sn + Kf, - R^)cn>lRH - (S^ - K„ - Ri^lc^’-OH (,,1^, 

r - /t+ ( 18a) 

= - V 08b) 

= Ry(Sj^ + Kj^) ■ (18c) 


and -H < z < 0. It should be noted that considerable caution is needed when the relation- 
ships are applied to field data, because radiance within canopies can vary drastically from 
place to place due to the random nature of the leaf spacing. The use of the solutions 
can be justified if the average radiances arc observed over large areas. In remote sensing 
applications of canopy reflectance models, reflected radiation is measured far above the 
canopy by airborne or satellite sensors, and thus more adequately represent the values of 
the canopy reflectance for the area viewed by the sensor. 
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REFLECTANCE/BIOMASS RELATIONSHIPS 
Tlie abilHy to monitor vegetation canopies using remote sensing techniques results from 
inherent or environmcn tally induced reflectance differences among plants and plant types 
during their growth cycle. In this section the reflectance/biomass relationships will be 
derived from the solution to the extended KM equations and their important properties 
will be discussed. 


The apparent diffuse reflectance formula for a vegetative canopy is obtained by letting 
z°Q from Eq. 16 as 

R « E^. (z«0)/Eg 

- Rb) Ry - tRy - Rk) ^ c (19a) 


R„ 


- R„ - (R„ - R„) 


(Sh t - R,,) (S^ - - (S^ - - R,,X5^ Kn>°' 

-f K, - R„ - (S^ - K^- Ri,) 


IK pH 


tl9b) 


g-2KpH ^ (r - RbRy) (Ry - R) 

(r - R Ry) (Ry - Rj,) 

_ (Sn + Kn - Ri,) (Sn - - R) 

(Sj^ + Kn - R) (Sn - - Rb) 


(19c) 

(19d) 


if Ry R|^ Tliis is the most general formula for the apparent canopy reflectance. It 
should be noted that the constant parameter D does not appear in this relationship: thus 
its expression is considerably less complex. Tliese exponential equations show asymptotic 
behaviors as the canopy height H (or biomass pH) approaches the two extreme values: 


and 


R- 

R- 


R., 


R,. 


as H *«- 00 (or pH 

as H — 0 (or pH 


oo) 

0 ). 


Tile asymptotic properties agree with most observations of canopy reflectance (Pearson, 
1973; Tucker, 1977; Deering, 1978). 
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Relationships Simplified 

It is higlily desirable to simplify complicated formulas for easier use and wider applicability, 
A couple of simplified relationships can be derived from Eq. 19. The general formula can 
be approximated by the simple exponential equation: 

R « Ry - (R, - Rjj) ^20) 

when r is sufficiently larger than one, since 0 < Ry, R, Rj^ < 1. Deviation of the apparent 
reflectance R by this equation from that by Eq, 19 is given approximately by 
(R^ - Ry)^ Ryexp(-2KpH)/r. Even if r is close to one, the approximation by Eq. 20 may 
be sufficient for practical use, since (Rj, - Ry)^RyCxp (-2KpH) < 1 in most cases. In field 
observations the uncertainty of refiectanee data is considerably large due to the presence 
of the direct solar radiation (Kriebel, 1976). Hence, such an approximation might bo 
justified in analysis of the biomass/refiectance data when bidirectional reflectance observa- 
tions are used as the diffuse reflectance data. In the equation the exponential term is 
the contribution of the background reflectance to the apparent canopy reflectance. It 
shows that the soil background effects will be significant when the difference between 
Ry and R^^ is large and biomass pH (or canopy height H) is small. 

Tlie asymptotic properties of the relationship are clearly seen at the two extreme values 
of the canopy biomass in Eq. 20. Tlie equation also indicates the ranges of the apparent 
canopy reflectance as 

Ry < R < Rb if 0< Ry < Rb < 1. 

and 

Ry > R > Rb if 1 > Ry ^ Rb > 0. 

Tliis working equation is simple in form and similar to many biological formulas possessing 
various asymptotic properties. For example, a similar relation has been employed to 
empirically fit spectroreflectanee and chlorophyll data (Pearson, 1973: Tucker, 1977). Such 
a simplicity would enhance the utility and acceptance of the relationship in agricultural 
applications. 
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Another commonly cited case is that and y_ = 7 ^. In this case Eq. 1 9c becomes 

e-2KpH = (1 - RuK)(K - R) , 

(1 - RRv) (Ry - Rb) 

since r=Y_/ 7 ^=l. A relationship equivalent to Eq. 21 has been derived by Allen and Rich* 
ardson (1968), who employed the leaf area index (LAI) instead of canopy biomass pH. Tlie 
relationship, wliich has been used to predict canopy LAI as a function of the apparent 
canopy reflectance at 0.4 pm, was depicted for shortgrass prairie vegetation (Bouteloua 
gracilis) and its asymptotic property svas also cited without elaboration by Oliver and 
Smith (1973). 

Characteristic Parameters 

llie values of R, Rjj, Ry and pH may be measured to some degree and, tlien the parameter 
K may be computed. However, a series of observations of apparent reflectance R and other 
canopy variables, such as H or pH, can lead to estimation of the other parameters: 

Ry and K or Kp. Once tlie values of the characteristic parameters are known for a growth 
cycle of vegetation, tlie plant bioma.ss can be assessed nondestructively by observing apparent 
canopy reflectances. It is possible to establish reflectance/biomass relationship curves for 
different crops or vegetation types by finding proper values of these characteristic parameters 
(Rjj, Ry and K or Kp, as well as r) for various agronomic and environmental conditions. 

Tlie value of r yields the first clue for the anisotropic property in canopy rencctance 
characteristics. 

Model Evaluation 

Tlie canopy reflectance relationships, Eqs. 19 tlirougli 21, were tested for the biomass/ 
reflectance data of 1) alfalfa and 2) shortgrass prairie vegetation. The alfalfa data were 
taken from seven experimental plots having different plant density. The different biomass 
levels were created througli selective tliinning within small plots (3m"), which were con- 
tained within a larger, uniform stand of alfalfa. Tlie canopy was approximately 45 - 50 cm 
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high, contained very little brown plant material, and was ready for the fall hay cutting. 

Alfalfa canopy spectral reflectances were acquired in a sampling mode (10 samples per 
small plot) using a two*channel, red and photographic infrared portable radiometer with 
spectral hands centered at about .68 and .80 urn, respectively. The observations were made 
under various sky and illumination conditions on three consecutive days (Table 1). 


Table 1. Reflectance of Alfalfa Canopies at Wavelength .68/im 
Observed on October 11-13, 1978 



OBSERVATION 



DRY BIOMASS (kg/ha) 



PLOT 

ILLUMINATION 

NOMINAL 

SUN 

ELEVATION 

l.D. 

DAY/TIME 

0 

70 

960 

1,650 

2,280 

3,660 

3,830 

CONDITION 

(deg.) 

A1 

11/ 9:13- 9:37 

.174 

.166 

.084 

.042 

.034 

.032 

,028 

Sunny 

32 

A2 

9:47.10:00 

.186 

.172 

.084 

.041 

.033 

,028 

.030 


36 

A3 

10:22.10:34 

.180 

.179 

.089 

.047 

.036 

.029 

,027 

H 

40 

A4 

11:03*11:19 

.189 

.169 

.091 

.030 

.033 

.028 

.027 

»l 

43 

A3 

12:00-12:14 

.186 

.172 

.m 

,047 

.031 

.027 

.027 

It 

44 

A6 

13:05.13:21 

.183 

.168 

.073 

.031 

.030 

.025 

.023 

H 

40 

A7 

14:01.14:ir 

.176 

.162 

,076 

.038 

.028 

,026 

.022 

II 

34 

no 

W:30.i3:07 

.173 

.139 

.070 

,039 

.033 

.023 

.026 

If 

27 

A9 

13:3(>-I3:43 

.189 

.164 

.073 

,067 

.030 

.026 

.024 

II * 

21 

B1 

12/13:30.14:24 

,178 

.160 

.073 

.047 

.032 

.023 

.026 

Hazy 

34 

B2 

14:41.14:36 

.183 

.160 

.078 

.039 

.033 

.030 

.029 

M 

28 

C 

13/12:31.13:03 

,233 

.212 

.108 

,081 

.047 

.033 

,035 

Hazy with 
high cirrus 

41 

D 

13/12:33.13:04 

.270 

.235 

.113 

.166 

,051 

.042 

.038 

Hazy with 
high cirrus 
and shaded 

41 


Under the clear sky conditions, measurements were made at nine different sample times 


during the day to examine the diurnal/sun elevation effects. Measurements were also made 
under more diffuse ligliting conditions of hazy skies and skies containing high thin cirrus 
clouds. One data set was also collected for a no-direct-sunliglit condition by artificially 
shading tire plant canopy. 


As the model was developed for diffuse illumination conditions, the reflectance data for 
the overcast day were potentially the most favorable for analysis of the model relationships. 
The clear sky conditions were considered to be a crucial test of whether the relationships 
were applicable to the anisotropic illumination, which is a deviation from the diffuse irra- 
diance assumption. 
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Tlie biomass/reflectance data of the sliortgrass prairie vegetation, dominated by blue grama 
grass (Bouteloua gracilis), were cnosen from the field data of 27 native grassland plots, as 
reported by Pearson and Miller (1972). Tlieir spectrorefifctances were measured, using a 
narrow band spectrometer of similar design as that used in the alfalfa field study. The cor- 
responding biomass data were tlrose of the total (green + dead) dry standing biomass clipped 
from each 1/4 square meter plot within the field-of-view of tlie spectrometer. 

Suitability of tlie model relationship to the refiectance/biomass data depends mostly on the 
parameter values. All the physics of canopy reflectance characteristics is depicted in terms 
of tliese bulk parameters of the model. Spectral reflectance of the background surface 
can be measured without difficulty if the areas of negligible canopy can be found in tlie 
study sites. Vegetation reflectance Ry miglit be approximated by extrapolation of data 
points when observations cover a wide range of canopy development stages or when observed 
data have “leveled off.” However, the parameter K (or the product Kp) is tlie resultant 
attenuation coefficient of tlie canopy as a half transparent medium, and its value can be 
obtained indirectly by evaluating Eq. 15, 16, 19 or 20 based on a series of observed data. 

For example, if all the otlier constant parameters (R^^ and Ry) are known, only one set 
of R and pH data will be enougli to determine K by Eq, 20. In the event that two or more 
available data sets may yield different values of K, the optimal solution for K is desired. 

Such an optimum estimate of au unknown parameter set can generally be found as a 
solution to the observation equations generated by a worxing equation for given data sets. 


In this investigation Rj^, Ry, and K were estimated such that their values yielded the best 

approximation of Eq. 19, 20 or 21 for given data set, that is, 

MinSW. jfj(R^, Ry, K;r; Bj, Rj)} 2, 
i 


where 


B- = (pH)j : Canopy biomass of the i-th obs. 


( 22 ) 
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( 23 ) 


! i , (r - RjRy) (Ry- Rfa) „-.2KB, 

(r - (Ry - 

1 _ (1 - RjRy) (Ry “ »~2KBj ^ (24) 

(1 - RjjRy) (Ry - Rj) 

1 -.fllfi! e-2KBi , (25) 

Ry - Rj 

Wj = weiglit for the i-tli obser\Mtion, 

and Rj is the observed spectral reflectance for the canopy of biomass B*. Eqs. 23, 24, and 
25 are equivalent to Eqs. 19, 21. and 20, respectively, if fj = 0. The best set of tlie param- 
eter estimates minimizes the weighted sum of the squared fj for the wliole data and was 
computed using tlie IMSL Subroutine ZXSSQ (IMSL Library, 1979). \Vlien the solution 
was searched by the subroutine, Eq. 23 was inefficient for computation and, hence, not 
used in the later analysis. 


Model and Observed Relationships 

The estimates of tlie characteristic parameters obtained by Eq. 24 or 25 (Table 2) produced 
curves clear •' depicting the observed relations between biomass and spectral reflectance 
(Figs. 2 and 3). No noticeable differences were found between the results from Eqs. 20 
and 21 (Table 2). For shortgrass prairie canopy data of Pearson and Miller no definitive 
best relationship could be drawn due to the large scatter of the data. In the optimum model 
solution tlie reflectance Ry of shortgrass prairie infinite canopy condition (i.e., that of the 
sufficiently large biomass as seen in the asymptotic character) was zero, wliich was out of 
the acceptable range for the solution, but it was certainly more realistic than the negative 
reflectance shown in the empirical Imear regression relationship (Fig. 3), which was achieved 
at only 4200 kg/ha. Tlie exponential formula, Eq. 20, worked favorably for analysis of the 
biomass/reflectance data — deviating by about 2.4% for biomass estimates and less than 0.2% 
for apparent reflectance estimates from the more complicated formula. 
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Table 2. Estimates of Reflectance Characteristic Parameters and Comparison of the Two 
Results Obtained by Eqs. 20 and 21 
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Fig, 2. Spectral Reflectance vs. Dry Biomass 
of Alfalfa Canopies at .68jum 


Fig. 3. Relationships of .68nm Reflectance 
vs. Total Dry Biomass for Short- 
grass Prairie Canopy Data of Pear- 
son and Miller (1972) 
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Tlie values of the canopy reflectance characteristic parameters, however, varied for differing 
illumination conditions, and so did the KM parameter (Table 2). Tlieir dependences on the 
sun angle seemed certain (Fig. 4). 
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Fig. 4. Reflectance Characteristic Parameters of Alfalfa Canopy at .68jum as a Function of 
the Solar Zenith Angle {B^) under Clear Sky Conditions 

Hence, bulk absorption and scattering coefficients of the vegetation canopy were a function 

of the insolation condition. Tlie coefficient of variation, of the ratio Ry/Rjj was 20 - 30% 

less than those of Ry and Rjj (Table 2), indicating that the ratios of the two reflectances, 

Ry and Rj^, remained fairly constant wliile the two individual parameters changed considerably 

depending upon illumination conditions. Although the illumination conditions were not ideal 

for the diffuse reflectance model, they are realistic remote sensing conditions, and tlie 

diffuse reflectance relationships were shown to be useful for accurate estimation of alfalfa 

and shortgrass prairie biomass utilizing measurements of plant canopy reflectance. 

SUMMARY 

Useful relationships were formulated to describe variations of the diffuse spectra) reflectance 
in terms of vegetation canopy variables, such as biomass. Tlie relationships were based on 
the solution of tire two-stream approximation of the radiative transfer equation. Out of 
tlie lengthy original expression of the diffuse reflectance formula, simple working equations 
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were derived by employing cliaracteristic parameters, wliicli are independ.ent of the canopy 
coverage and identifiable by field observations. The typical asymptotic nature of reflectance 
data that is usually observed in biomass studies was clearly explained. It also established the 
range of expected apparent canopy reflectance values. 

A procedure to estimate reflectance cliaracteristic parameters was described for practical 
applications of the relationships. The simplified exponential formulas accurately depicted 
the obsei^'ed relationships between biomass and spectral rcfiectance. They were shown to 
be useful for accurate estimation of alfalfa and shortgrass prairie biomass utilizing measure- 
ments of plant canopy reflectance. 


15 



REFERENCES 


Allen, W.A., and AJ. Richardson, 1968, Interaction of Light with a Plant Canopy, 

J. Opt. Soc. Am,, Vol. 58, No. 8, pp. 1023-1028. 

Breece, H.T. (Ill) and R.A. Holmes, 1971, Bidirectional Scattering Characteristics of Heaithy 
Green Soybean and Corn Leaves in YiyQ,Appl Opt., Vol. 10, No. 1, pp. 119-127. 

Deering, D.W. 1978, Rangeland Reflectance Cliaracteristlcs Measured by Aircraft and Space- 
craft Sensors, Ph.D. diss., Texas A&M University, College Station, Te,xas, 316 p. 

Duggin, M.J. 1977, Likely Effects of Solar Elevation on the Quantification of Changes in 
Vegetation with Maturity Using Sequentlai Landsat Imagery, App/. Opt,, Vol. 16, No. 3, 
pp. 521-523. 

Gates, D.M. 1965, Characteristics of Soil and Vegetated Surfaces to Reflected Emitted Radia- 
tion, pp. 573-600, Jn: Proc, Tliird Int, Symp. on Remote Sens, of Environ., Univ. Michigan, 
Ann Arbor. 

Holben, B.N., C.J. Tucker, and C.J. Fan, 1980, Spectral Assessment of Soybean Leaf Area 
and Leaf Biomass, Photogram, Engr, and Remote Sens. Vol. 46, No. 5, pp. 651-656. 

IMSL Library, Reference Manual, 1979, IMSL, Inc., Houston, Texas. 

Kimes, D.S., J.A. Smith, and K.J. Ranson, 1980, Interpreting Vegetation Reflectance 
Measurements as a Function of Solar Zenith Angle, Photogram, Engr, and Remote Sens. 

(in press). 

Kortiim, G. 1969, Reflectance Spectroscopy, Springer-Verlag New York, Inc., 366 p. 

Kriebel, K.T. 1976, On the Variability of the Reflected Radiation Field Due to Differing 
Distributions of the Irradiation, Remote Sens, of Environ., Vol. 4, pp. 257-264. 


16 



Kriebel, K.T. 1978, Measured Spectral Bidirectional Reflection Properties of Four Vegetated 
Surfaces, Appl. Opt., Vol. 17, No. 2. pp. 253-259. 

Loomis, W.E. 1965, Absorption of Radiant Energy by Leaves, £co/., Vol. 46, pp. 14-17. 

Oliver, R.E., and J.A. Smith, 1973, Vegetation Canopy Reflectance Models, Final Report to 
U.S, Army Research Office - Durham, DA-ARD-D3 1-1 24-7 1-G 165, Colorado State Univ., 
Fort Collins, 82 p. 

Park, J.K. 1980, A Soil Moisture Reflectance Model in Visible and Near IR Bands, presented 
at Int. Symp. of Machine Proc. of Remotely Sensed Data and Soil Infor. Syst. and Remote 
Sens, and Soil Survey, June 3-6, Purdue Univ., West Lafayette (NASA TM 80701). 

Pearson, R.L. 1973, Remote Multispectral Sensing of Biomass. Ph.D. diss. Colorado State 
Univ,, Fort Collins, 180 p. 

Pearson, R.L., and L.D, Miller, 1972, Remote Mapping of Standing Crop Biomass for Esti- 
mation of the Productivity of the Shortgrass Prairie, Pawnee National Grasslands, Colorado, 
pp, 1355-1379, In; Proc. Elglith Int. Symp. on Remote Sens, of Environ., Univ. Michigan, 
Ann Arbor. 

Thompson, D.R., and O.A. Wehmanen, 1979, Using Landsat Digital Data to Detect Moisture 
Stress, Photogmm. Engr, and Remote Sens., Vol. 46, No. 5, pp. 651-656. 

Tucker, C.J. 1 977, Asymptotic Nature of Grass Canopy Spectral Reflectance, ^ippl. Opt., 
Vol. 16, No. 5, pp. 1151-1156. 

Tucker, C.J. 1979, Red and Photographic Infrared Linear Combinations for Monitoring 
Vegetation, Remote Sens, of Environ,, Vol. 8, No. 2, pp. 127-150. 


17 



Tucker, CJ. 1980, A Critical Comparison of Remote Sensing and Other Nondestructive 
Biomass Estimation Metliods, Forage ScL (in press) (NASA TM 80607). 

Wendlandt, W.W., and H.G. Heelit, 1966, Reflectance Spectroscopy, Interscience Publisher, 
298 p. 

VViegand, C.L, A.J. Richardson, and E.T. Kanemasu, 1979, Leaf Area Index Estimates for 
Wheat from Landsat and their Implications for Evapotranspiration and Crop Modeling, 
Agron. Vol. 71, pp. 336-342. 


18 



RELATIONSHIPS BETWEEN DIFFUSE REFLECTANCE AND VEGETATION CANOPY 


VARIABLES BASED ON THE RADIATIVE TRANSFER THEORY* 

John K. Park ami Donald W. Deerino 
Earth Survey Applications Division 
NASA/Goddard Space Flight Center 
Greenbelt, MD 20771 

BIOGRAPHICAL SKETCH 

John K, Park received his Ph,D. in Civil Engineering from Colorado State University, Fort 
Collins, Colorado, in 1979, specializing in cluster analysis of multispectral imagery data. 

His interests are in pattern recognition and mathematical modeling of soil moisture and 
vegetation canopy reflectance for assessment of crop water demands by remote sensing 
techniques. He is currently a postdoctoral Resident Research Associate in the Hydrological 
Sciences Branch of Goddard Space Flight Center, Greenbelt, Maryland. 

Donald W. Deering is a physical scientist in the Earth Resources Branch of Goddard Space Flight 
Center, Greenbelt, Maryland, where he specializes in fundamental field and aircraft studies di- 
rected toward future satellite sensor definition for the remote sensing of vegetation. He received 
his M.S, in 1972 from Texas Tech University and his Ph.D. in 1978 from Texas A&M University. 
Botli degrees were in Range Science. From 1972 to 1978 he was a Research Associate with 
the Remote Sensing Center at Texas A&M University conducting research in rangeland remote 
sensing. His research at Goddard Space Fllglit Center has been concentrated on bidirectional 
reflectance and sun elevation effects on sensor response and vegetation parameter assessment. 


* This paper Is to be presented at the 1881 ASP-ACSM Annual Convention, Washington, D.C„ February 23-27, 1981. 


19 



